The Adam-Gibbs relation for glass-forming liquids in 2, 3 and 4 dimensions 
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The Adam-Gibbs relation between relaxation times and the configurational entropy has been 
tested extensively for glass formers using experimental data and computer simulation results. Al- 
though the form of the relation contains no dependence on the spatial dimensionality in the original 
formulation, subsequent derivations of the Adam-Gibbs relation allow for such a possibility. We 
test the Adam-Gibbs relation in 2, 3, and 4 spatial dimensions using computer simulations of model 
glass formers. We find that the relation is valid in 3 and 4 dimensions. But in 2 dimensions, the 
relation does not hold, and interestingly, no single alternate relation describes the results for the 
different model systems we study. 

A central theme in the study of glass forming liquids is a satisfactory understanding of the behaviour of relaxation 
times as the glass transition is approached, showing temperature dependence typically stronger than the Arrhenius 
law: t(T) = t(oo) exp[-^^] observed at high temperatures GH3. The Adam-Gibbs (AG) relation which is of 
central importance in glass forming liquids, explains the behaviour of the relaxation time i.e. dynamics in terms of the 
configurational entropy i.e. thermodynamics. It also forms the basis of more sophisticated theories of glass transition 
connecting dynamics to thermodynamics e.g. the random first order transition theory (RFOT) In these theories 

the spatial dimensionality (D) appears explicitly in the relationship between dynamics and thermodynamics. However, 
the large number of experimental and numerical studies wherein the AG relation has been tested 1Q- 17 1 have all been 



in 3 dimensions. The present study aims to critically examine whether the AG relation is valid in different spatial 
dimensions or gets generalized in a D dependent way. 

The AG relation is based on the picture that relaxation in glass forming liquids occurs through the collective rear- 
rangement of "cooperatively rearranging regions" (CRR). The CRRs define a minimum size of groups of rearranging 
particles (atoms, molecules etc. depending on the nature of the glass former) such that smaller groups of particles are 
incapable of rearrangement independently of their surroundings. Adam and Gibbs argued that the configurational 
entropy (the entropy associated with the multiplicity of distinct arrangements of particles, obtained by subtracting 
a "vibrational" component from the total entropy) per particle S C (T) of a liquid varies inversely as the size of the 
CRR, z(T), since the configurational entropy per CRR, S*, is roughly independent of temperature: S C (T) = -Jjy 
The further assumption that the free energy barrier for a rearrangement is proportional to the size of the CRR 
(AG = zSfi, 5jjL = chemical potential barrier per particle) results in the Adam-Gibbs relation: 



t(T) = r(oo)exp 
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The above relation is obtained independent of reference to the spatial dimensionality of the system and is hence 
expected to be same in all spatial dimensions. 

A rationalization of the AG relation, based on more detailed considerations of possible activated relaxation mech- 
anisms, is offered by the RFOT, which has recently been discussed by many authors specifically in the context of 



growing length scales associated with the glass transition [1 81-42-4 1 



In the "mosaic" picture @,0], the liquid is divided into metastable regions of characteristic size £(T). The transition 
from one metastable state to others is hindered by the cost of surface free energy ( AG S oc Y^ e , where 9 < D, where 
D is the spatial dimension and Y is the surface tension) and driven by the possibility of sampling an exponentially 

large number of other metastable minima (free energy gain AG& oc TS C ^ D ). There is a characteristic cross-over length 

i 

£*(T) oc ^f^r) above which a liquid samples all metastable states and below which the liquid is trapped in one 
of the metastable states. The characteristic length £(T)* diverges as Sc(T) —> 0. Assuming in general that the free 
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energy barrier to relaxation varies (following the notation of @) as AG(T) oc £(T)^, one obtains 



t(T) = r(oo) exp 



.4 



l(TS c )- 



(2) 



where A has weak T dependence, and ip = 9 if AG(T) is calculated as the free energy barrier obtained from the 
surface and bulk contributions above (AG = AG& + AG S ). The mosaic picture thus contains an explicit dependence 
on the spatial dimesion, and one recovers the Adam Gibbs relation only when jjzrg — 1- The original AG formulation 
is equivalent to assuming ip = D, 9 = 0, whereas it was argued in Q that 9 = ip = D/2, with the AG relation as the 
result in both cases. In the latter case a dimension independent AG relation is predicted in spite of arguments 
that treat spatial dimensions explicitly. 

However, attempts to estimate the exponents either numerically or from experimental data are not conclusive. By 
direct simulation in a model liquid (D = 3), Cammarota et al. obtained 6 = 2, ip = 1 [22| which, although consistent 
with the AG relation, differs from both the AG and RFOT values. The simulation study by Karmakar et al. 2^. 241 
obtained 9 ~ 2.3 and hence ip ~ 0.7 using the condition 



D- 



= 1. An extensive study by Capaccioli et al. [20j of 
experimental data for 45 glass-forming liquids show a lot of variation in both 9 and ip, with best fit estimates in the 
range of 9 ~ 2 — 2.15, and ip ~ 0.85 — 1. These two studies are consistent in suggesting values of 9 > 2, ip < 1, but 
both use measures of dynamical heterogeneity to extract a length scale. Since a priori one must expect the mosaic 
length to be distinct from the heterogeneity length scale, the implications of these estimates arc not clear. 

In view of the considerations above, the possibility of an explicit dependence of the AG relation on the spatial 
dimension (and within the RFOT framework, the values of the exponents 9 and ip) merits investigation. In this letter, 
we address the question of the dependence on spatial dimension of the AG relation by studying model liquids in 2, 
3 and 4 dimensions, and evaluating the relationship between relaxation times and the configurational entropy. The 
evaluation of the RFOT exponents 9 and ip, while relevant to our investigation, requires the calculation of the mosaic 
lengthscale. Ways of obtaining lengthscales relevant to glassy dynamics that have been discussed in the literature (see 



e -g-j 21, HI IS 43 1) tend to be computationally very demanding 42 1. We have therefore not attempted to evaluate 



the RFOT exponents in this work. 

One of the well-known systems in which the AG relation is shown to be valid [H , 23 , 2|| is Kob- Andersen model 
(KA) [26j . which is an 80 : 20 binary mixture of particles interacting with a Lennard- Jones (LJ) potential, with LJ 
parameters cab / ^aa = 1-5, cbb/caa = 0.5, oab/o~aa — 0.80, ubb/o'aa — 0.88. We truncate the interaction potential 

at 2.5cr a ^ (details as in 15|). Units of length, energy and time scales are oaa^aa and J aA ^ mAA respectively. We 



have studied the KA model in two (2D), three (3D) and four (4D) spatial dimensions. In addition in two dimensions 
we have studied the following liquids: (1) KA model at 65:35 composition [StJ which is denoted as the modified KA 



(MKA) model and (2) 50:50 binary mixture of repulsive soft spheres (i?10 model, V(r) 



-10 



) (see 



for details of 

the model. Like the KA models, we use reduced units defined in terms of the energy scale of the model and the size 
of the large particles). We have performed molecular dynamics simulations in the canonical (NVT) ensemble, using 
the constant temperature algorithm of Brown and Clarke [29j . Simulations were done at a fixed number density p 
(p = 1.2 for the KA model in two and three dimensions, p = 1.6 in four dimensions, p = 1.2 for the MKA model 
and p = 0.85 for the R10 model). Integration time steps were in the range £ [0.001,0.006] depending on T. Run 
lengths at each T were in excess of 100r Q (the relaxation time r Q is defined below) (e.g., at the lowest T in 4D the 
total runlcngth is ~ 6 x 10 8 MD steps). The range of relaxation times accessed is upto O(10 4 ) in three and four 
dimensions and upto O(10 5 ) in two dimensions. Typically, 3 (1 — 3 in 4D, 3 — 5 in 3D and 2D) independent runs are 
performed at each T. In order to assess the influence of system size, we show results for different system sizes in 2 
and 3 dimensions (indicated in appropriate places) although we do not discuss system size effects in detail here. 

As a measure of dynamics, we have studied a two- poin t time correlation function, the overlap function q(t) and its 
fluctuation related to dynamic susceptibility (X4(*)) |23l 30l- 32 1 : 



<q(t)> = < drp(r,t )p(r,t + t ) > 



N 



< J2w(\n(t ) - n(t +t)\) > 



X4(t) = ^(q(t) 2 )-(q(t)) 2 
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where p(r,t ) is the space and time dependent particle density, w(r) is the window function: w(r) = l,r < 
a and zero otherwise (a — 0.3 in 3D) and averages are obtained over initial times to as well as independent sam- 
ples. Relaxation times are estimated both from the condition q(r)/N = 1/e (r Q ) and the characteristic time where 
Xi(t) is maximum (T4). As we find these two times to be proportional to each other [23|, we report only r a . 

The configurational entropy (S c ) per particle, the measure of the number of distinct local energy minima, is 
calculated [33j by subtracting from the total entropy of the system the "vibrational" component: 

S C (T) = StotaiiT) — S V ib(T) 

The total entropy of the liquid is obtained via thermodynamic integration [33| . The vibrational entropy is calculated 
by making a harmonic approximation to the potential energy about local energy minima (termed "inherent strcture"s) 
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We evaluate the configurational entropy below the onset temperature 
Arrhenius to non-Arrhenius behaviour of the relaxation time occurs TBI 
in the 4D KA model. Fig. 1(b) shows the temperature variation of the inherent structure energy which exhibits a 

ET 



Tonset) across which a cross-over from 
2l. In Fig. ~ 



1(a) we show this cross-over 



1 /T temperature dependence below T onset [15( . In Fig. 1 (c) we show the T dependence of the configurational entropy 
for the 4D KA model. As T is lowered TS C goes to zero linearly with T. We estimate the Kauzmann temperature 
(Tr-) at which the extrapolated configurational entropy vanishes to be Tk = 0.53. 

In Fig. [5] we show the AG plot for the KA model in different spatial dimensions. We see that in 4D, the AG relation 
is valid over several orders of magnitude of change in r a . As reported in earlier works, the AG relation is also valid in 
3D, and the data presented in Fig. 2(b) indicates a mild system size dependence in the slopes C (Eq. [1]). However, 
in two dimensions [Fig: 



2(c) 



we see deviations from the AG relation for low temperatures and small system sizes. 
It has been reported recently that the 2D KA model is prone to orientational ordering 27j, which we confirm (data 
not shown). It is therefore not clear to what extent these observed deviations are significant. In order to asertain 
the dependence of r Q on S c without the presence of orientational ordering, we study two other previously studied 2D 
models which do not exhibit significant ordering, namely the KA model with composition 65 : 35 (MKA) j27j and 
the 50 : 50 binary mixture of repulsive soft spheres (RIO model, V(r) ~ r~ 10 ) [28|. Figs 3(a) and |3(b)j show that for 
these systems, r a does not obey the AG relation. 

We note that the nature of deviation is different in these two cases. Fitting r a to the generalized AG form 
t = t cxp ((t§-) q ) , we find a = 0.43 for N = 500, a = 0.60 for N = 2000 and a = 0.67 for N = 10000 (i. e. a < 1) 

for the MKA model, whereas for the i?10 model we find a = 2.1 (i. e. a > 1). 

The MKA model has both attractive and repulsive interactions, whereas the R10 model has purely repulsive 



interactions. It has recently been suggested [39j that the exponent ip in the relation r(T) = r(oo) exp (see 
discussion preceding Eq. ^ depends on the nature of the interactions. This possibility suggests an explanation for the 
results we observe which is under investigation. However, we note that in three dimensions, the AG relation has been 
verified in both attractive models (e.g. Kob- Andersen [l5|, Lewis- Wahnstrom OTP [l6| and Dzugutov liquid [4(|) 



and repulsive models (e.g. repulsive soft spheres |41(). We also note that the exponent a above gets larger for larger 
system sizes for the MKA model. Although N = 10000 is a sufficiently large system as far as any relevant length 
scale in this model is concerned, finite-size effects as the orgin of the deviation from the AG relation in these systems 
cannot be completely ruled out at present. Finally we note that the behaviour described here remains qualitatively 
the same if we use inverse diffusivities instead of a-relaxation times. Such a comparison will be presented elsewhere 
[iij in the context of the breakdown of the Stokes-Einstein relation. 

In summary, we examine the validity of the AG relation in different spatial dimensions by studying model liquids 
via computer simulations, for the first time in 2 and 4 dimensions. The Adam-Gibbs relation is valid in four and 
three dimensions, but is not obeyed in 2 dimensions and the nature of the deviation from the AG relation depends 
on the details of the interaction between particles. This is an unexpected result. Although an understanding of this 
observation is lacking at present, they present constraints that a successful theoretical explanation of slow relaxation 
must meet. Such a lack of universality is unexpected and surprising since both Adam-Gibbs and RFOT theories 
describe the glass transition as a thermodynamic phase transition. 
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FIG. 1: 1(a) Temperature dependence of the relaxation time r a of the 4D KA model displaying a cross-over from Arrhenius to 
non-Arrhenius behaviour at T onset ~ 1. The solid line is an Arrhenius fit above T onse t. |l(b) shows the average inherent structure 
energy vs. T displaying a 1/T dependence below T onset . 1(c) Kauzmann temperature for the 4D KA model, Tk = 0.53, is 
obtained from the condition TkS c (Tk) = 0. 



[2] C. A. Angell, K. L. Ngai, G. B. McKenna, P. F. McMillan, S. W. Martin, J. App. Phys., 88, 3113 (2000). 
[3] P. G. Debenedetti and F. H. Stillinger, Nature 410, 259 (2001). 
[4] G. Adam and J. H. Gibbs, J. Chem. Phys. 43, 139 (1965). 

[5] T. R. Kirkpatrick, D. Thirumalai, and P. G. Wolynes, Phys. Rev. A 40, 1045 (1989). 

[6] G. Parisi and M. Mezard, J. Chem. Phys. Ill (1999). 

[7] X. Xia and P. G. Wolynes, Proc. Natl. Acad. Sci. 97, 2990 (2000) 

[8] V. Lubchenko and P. G. Wolynes, Annu. Rev. Phys. Chem. 58, 235 (2007); 

[9] J. Bouchaud and G. Biroli, J. Chem. Phys. 121, 7347 (2004). 
[10] R. J. Greet and D. Turnbull, J. Chem. Phys. 47, 2185 (1967). 
[11] R. Richert and C. A. Angell, J. Chem. Phys. 108, 9016 (1998). 
[12] R. J. Speedy, Mol.Phys. 95, 169 (1998). 
[13] G. P. Johari, J. Chem. Phys. 112, 8958 (2000). 

[14] A. Scala, F. W. Starr, E. La Nave, F. Sciortino and H. E. Stanley, Nature (London) 406, 166 (2000). 
[15] S. Sastry, Nature 409, 164 (2001). 

[16] S. Mossa, E. La Nave, H. E. Stanley, C. Donati, F. Sciortino, and P. Tartaglia, Phys. Rev. E 65, 041205 (2002). 
[17] C. M. Roland, S. Capaccioli, M. Lucchesi and R. Casalini, J. Chem. Phys. 120, 10640 (2004). 



5 




(TS )"' / E 

1 C' ; 

(c) 



FIG. 2: Adam-Gibbs plots (r a "vs. (TSc) -1 ) for the KA model in 4, 3, and 2 dimensions. The AG relation holds in three 
and four dimensions but deviations are observed in two dimensions at low T and small system sizes. The system size in 4D is 
N = 1500. Lines are fits to the AG relation (Eq. [TJ (for the largest system size in the case of 2D). 
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FIG. 3: Adam-Gibbs plot for the MKA and RIO models. Systematic deviations from the AG relation are seen for both models, 
but with opposite deviations from linearity. For the MKA model, the lines are fits to the generalized AG relation of the form 

r a = r exp ((^§7)") (see also Eq. 0[a = 0.43 for N = 500 (dashed), a = 0.60 for N = 2000 (dotted line) and a = 0.67 

for N — 10000 (solid line)]. For the iilO model, the dashed line is a fit to the AG relation and the solid line is a fit to the 
generalized AG relation with a — 2.1. The system size is N — 2048 for the i?10 model. 
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